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1 Introduction 

In this paper we study the classical consumption-investment model with infinite 
horizon in the presence of transaction costs. Our aim is to extend the results of [18] 
to the case where the price evolution is given by a geometric Levy process. Namely, 
we show that the Bellman function is a viscosity solution of the corresponding 
Hamilton-Jacobi-Bellman equation. We also prove a uniqueness theorem for the 
latter. 

Mathematically, the consumption-investment problem with transaction costs 
is a regular-singular control problem for a linear stochastic equation in a cone. 
Its specificity is that the Bellman function may not be smooth and, therefore, 
one cannot use the verification theorem (at least, in its traditional form) because 
the Ito formula cannot be applied. Nevertheless, one can show that the Bellman 
function is a solution of the HJB equation in viscosity sense. Though the general 
line of arguments is common, one needs to re-examine each step of the proof. In 
particular, for the considered jump-diffusion model, the HJB equation contains an 
integro-differential operator and the test functions involved in the definition of the 
viscosity solution must be “globally” defined. It seems that already in 1986 H.M. 
Soner noticed that the control problems with jump parts can be considered in the 
framework of the theory of viscosity solutions, [26] . m- 

There is a growing literature on extension of the concept of viscosity solutions 
to equations with integro-differential operators, see, e.g., [23], [2], [22], 0 , 0 , 0 , 
[3]. There are several variants of the definition of viscosity solution. Our choice is 
intended to serve the model with a positive utility function. The definition can be 
viewed as a simplified version of that adopted in m- 

A rather detailed study of consumption-investment problems under transac¬ 
tion costs when the prices follow exponential Levy processes and the investor is 
constrained to keep long positions in all assets, money included, was undertaken 
in papers by Benth et al. EH and m■ Our geometric approach seems to be more 
general than that of the mentioned papers where the authors consider a ’’paramet¬ 
ric” version of the stock market, with transactions always passing through money 
(i.e. either “buy stock” or “sell stock”). A more important difference is that in our 
setting the investor may take short positions as was always assumed in the classical 
papers ed, m , [25] . If short positions are admitted, the ruin may happen due to 
a jump of the price process. That is why the natural, ’’classical”, setting considered 
here leads to a different HJB equation of a more complicated structure. Follow¬ 
ing the ideas from the paper [TS] we derive the Dynamic Programming Principle 
splitted into two separate assertions. Though it is the principal tool which allows 
to check that the Bellman function is a viscosity solution of the LIJB equation, it 
is rarely discussed in the literature (and even taken for granted, see, e.g., in 0, 

m , m)- 

The main results of the paper is Theorem 110.11 claiming that if the Bellman 
function is continuous up to the boundary then it is a viscosity solution of the 
HJB equation and the uniqueness theorem for the Dirichlet problem arising in 
the model, Th. 111.21 We formulate the latter in terms of the Lyapunov function, 
an object that is defined in terms of the truncated operator, in which the utility 
function is not involved. Its introduction allows us to disconnect the problems of 
the uniqueness of a solution and the existence of a classical supersolution. 
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Probably, the most important result of the paper is the uniqueness theorem for 
the HJB equation with a non-local operator. In contrast to the methods developed 
in [?] which are based on (very technical) extensions of the Ishii lemma we use the 
latter in its original and (very transparent) formulation. 

The structure of the problem is the following. In Sections 2 and 3 we introduce 
the model dynamics and describe the goal functional providing comments on the 
concavity of the Bellman function W. In Section 4 we show that the Bellman 
function, if finite, is continuous in the interior of the solvency cone. In Section 
5 we give a formal description of the HJB equation. Sections 6 and 7 contain a 
short account of basic facts on viscosity solutions for integro-differential operators. 
In Section 8 we explain the role of classical supersolutions to the HJB equations. 
Section 9 is devoted to the Dynamic Programming Principle. In Section 10 we use 
it to show that the Bellman function is the solution of our HJB equation. Section 
11 contains a uniqueness theorem formulated in terms of a Lyapunov function. In 
Section 12 we provide examples of Lyapunov functions and classical supersolutions. 


2 The Model 

Our setting is more general than that of the standard model of financial market 
under constant proportional transaction costs. In particular, the cone K is not 
supposed to be polyhedral. We assume that the asset prices are geometric Levy 
processes. Our framework appeals to a theory of viscosity solutions for non-local 
integro-differential operators. 

Let Y = (Y t ) be an Revalued semimartingale on a stochastic basis (12, T, F, P ) 
with the trivial initial cr-algebra. Let K and C be proper closed cones in R d such 
that C C int K 0. Define the set A of controls 7r = (B, C ) as the set of predictable 
Revalued cadlag processes of bounded variation such that, up to an evanescent 
set, 

B € —K, c € C. (1) 

Here B denotes a (measurable version) of the Radon-Nikodym derivative of B 
with respect to the total variation process ||B||. The notation C has a similar 
sense. Though models with arbitrary C is of interest, we restrict ourselves in the 
present paper by considering consumption processes admitting intensity. To this 
end, we define A a as the set of controls n with absolutely continuous components 
C such that the increment Co = 0. For the elements of A a we have c := dC/dt € C. 
The controlled process V = V x ’ n is the solution of the linear system 

dVj^vtdYj + dBt-dCi, V 0 i -=x i , i = l,...,d. (2) 

In general, AVq = ABq is not equal to zero: the investor may revise the portfolio 
when entering the market at time zero. 

The solution of m can be expressed explicitly using the Doleans-Dade expo¬ 
nentials 

£ t (y J ) = e W-(i/2)<T“=)t JJ(i + AYs)e~ AY * . (3) 

S<t 

Namely, 

V^ = St(Y i )x i + et(Y i ) [ S-'iY^idBi-dCi), i = l,...,d. (4) 

J[o,t] 
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We introduce the stopping time 

8 = e x ’ v := inf{t : V?'* $ int K}. 


For a; £ int K we consider the subsets A x and AS of “admissible” controls for which 
7r = /[Q gi, T ]7r and {V- + AB £ int A'} = {V_ £ int A}. In financial context, 8 is 
the time of ruin. When V x ’ n leaves the interior of the solvency cone the control 
of the portfolio and the consumption stops. The process V given by 0 , continues 
to evolve after the time 8 but for us the relevance has only the stopped process 

yX.TT.O 

It is natural to assume that the process V does not leave the interior of K due 
to a jump of B: the investor is reasonable enough not to ruin himself by making a 
too expensive portfolio revision. 

The important hypothesis that the cone K is proper, i.e. K n ( —K ) = {0}, 
or equivalently, int K* ^ 0, corresponds to the model of financial market with 
efficient friction. In a financial context K (usually containing R+) is interpreted as 
the solvency region and C = ( Ct ) as the consumption process; the process B = ( Bt) 
describes accumulated fund transfers. In the “standard” model with proportional 
transaction costs (sometimes referred to as the model of currency market) 

K = cone{(l + A u )ej — ej, e*, 1 < i, j < d} 

where \ lJ > 0 are transaction costs coefficients, see Section 3.1 in the book [T9] 
for details and other examples. 

The process Y represents the relative price movements. If S 1 is the price process 
of the ith asset, then dS\ = S\_dYf and S\ = S l 0 St{Y' 1 ). Without loss of generality 
we assume that Sq = 1 for all i. In this case Y l is the so-called stochastic logarithm 
of S l . The formula 0 can be re-written as follows: 


Vf = S\x 1 + si 


l 


[o,t] s l- 


( dBl-d & s ), i = l,..., d. 


We shall work assuming that 


Y t = ffi + Sw t + 


/;/ z(p(dz, dt) — q(dz, dt )) 


( 5 ) 

( 6 ) 


where p £ R^, w is a m-dimensional standard Wiener process and p(dz,dt) is 
a Poisson random measure with the compensator q(dz,dt) = n(dz)dt such that 
FI ( dz ) is a measure concentrated on ] — 1, oo[ d . Note that the latter property of the 
Levy measure corresponds to the financially meaningful case where S l > 0. For 
the m x d-dimensional matrix S we put A = SS*. We assume that 


J ( \z\ 2 A \z\)n(dz) < oo 

and this assumption validates the formula ©: by definition, 


( 7 ) 
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where the first integral is defined as a stochastic one, while the second and the 
third are the usual Lebesque integrals, both finite (a.s.). 


Notation. By typographical reasons we shall use the notation D x instead of 
common diaga: for the diagonal operator (or matrix) generated by the vector 
x = (x 1 ,..., x d ), i.e. 


r-x /Id d d\ 

D X Z = (X Z Z ). 


The system (J2J) can be written in the integral vector form as follows: 


V t = x + 


[ D Vs _(nds + EdW s )+ I 
Jo Jo 


Dv a - Z (p(dz,ds) - q(dz,ds)) + B t -Ct- (8) 


It is important to note that the jumps of Y and B cannot occur simultaneously. 
More precisely, the process |ziB||Z\F| is indistinguishable of zero. Indeed, for any 
e > 0 we have, using the predictability of the process AB = B — B -, that 


e^2\ab s \\ay s \i {1 

AY a \>e} 

s>0 



\ AB s\I{\z\>e}\ z \p(dz,ds) 


= E 



\ AB AW{\z\>e} n {dz)ds = 0 


because for each uj the set (s : AB s (ui) ^ 0} is at most countable and its Lebesgue 
measure is equal to zero. Thus, the process \AB\\AY\In^y\>e} i s indistinguishable 
of zero and so is the process |AB||zAY|. 

It follows that ABg = 0. Since the predictable process Iiv_edK}I[0,0] has at 
most countable number of jumps, the same reasoning as above leads to the con¬ 
clusion that I{v_^dK}\^Y\I[o,9] is indistinguishable of zero. This means that 9 is 
the first moment when either V or V- leaves int K. This property will be used in 
the proof that W is lower semicontinuous on int K. 

In our proof of the Dynamic Programming Principle (needed to derive the HJB 
equation) we shall assume that the stochastic basis is a canonical one, that is the 
space of cadlag functions and P is a measure under which the coordinate mapping 
is the Levy process. 


3 Goal Functionals and Concavity of the Bellman Function 

Let JJ : C —> R+ be a concave function such that 1/(0) = 0 and U(x)/\x\ —» 0 as 
\x\ —» oo. With every 7r = ( B,C ) £ we associate the “utility process” 

r-tf\9 

J f := / e~d s U{c s ) ds , t > 0, 

Jo 

where /3 > 0. We consider the infinite horizon maximization problem with the goal 
functional EJ and define its Bellman function W by 

W(x):= sup LJj,, zeintAL (9) 

irG A% 

Since Aa 1 C Aa 2 when X2 — x\ £ K, the function W is increasing with respect 
to the partial ordering >x generated by the cone K. 
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If TVi , i = 1,2, are admissible strategies for the initial points Xi , then the strategy 
A 7 ti + (1—A)7 T2 is an admissible strategy for the initial point A®i+(1— X)x 2 , A £ [0,1], 
laying on the interval connecting x\ and x -2 ■ In the case where the relative price 
process Y is continuous, the corresponding ruin time for the process 



( 10 ) 


dominates the maximum of the ruin times for processes V Xi,ni . The concavity of 
u implies that 



( 11 ) 


and, hence, the function W is concave on int K. 

Unfortunately, in our main case of interest, where Y has jumps, the ruin times 
are not related in such a simple way since the short positions are allowed. It is easy 
to give examples of trajectories such that 9 Xl ’ ni = 9 Xx i+P- x ) x 2 , x ni+(i-\)n 2 < ^ 
while 9 X2 ’ n2 = oo and the relations m and ED do not hold. Therefore, we cannot 
guarantee, by the above argument, that the Bellman function is concave. Of course, 
these considerations show only that the concavity of W cannot be obtained in a 
straightforward way as for a model based on continuous price process but it is not 
excluded. 

The concavity of the Bellman function W is not a property just interesting per 
se. The classical definition of viscosity solution, as was given by the famous “User’s 
guide” m, requires the continuity of W. On the other hand, a concave function 
is continuous in the interior of its domain (and even locally Lipschitz), see, e.g., 
[5]. Of course, the model must contains a provision which ensures that W is finite. 
But the latter property, in the case of continuous price processes implies, that W 
is continuous on int K. In the case of processes with jumps one needs to analyze 
the continuity of W using other arguments. 

In the next section we show that the finiteness of W still guarantees its conti¬ 
nuity in the interior of K. We do this using the following assertion. 

Lemma 3.1 Suppose that W is a finite function. Let x £ int A'. Then the function 
A W(Ax) is right-continuous on R+. 

Proof. Let A > 0. Then An £ Aa X if and only if n £ AS- For a concave function U 
with 1/(0) = 0 we have, for any e > 0 the inequality U(c) > (1 + e) _1 [/((l + e)c). 
Hence, for an arbitrary strategy n £ A% we have, for 9 = 9 X,V = 5 »( 1 + e ) x >( 1 + e ) 7r ) 
that 



It follows that W((l + e)x) < (1 + e)W{x). Since W(x) < W((l + e)x), we infer 
from here that A IU(Ax) is right-continuous at the point A = 1. Replacing x by 
Ax we obtain the claim. □ 

If U is a homogeneous function of order 7 with 7 £]0,1[, i.e. U{Ax) = A 1 U{x) 
for all A > 0, * £ K, then W(A*) = A' r W(x). Thus, the function A W{Ax) is 
concave and, therefore, continuous if finite. 
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Remark 1 . In financial models usually C = R+ei and <r° = 0, i.e. the only first 
(non-risky) asset is consumed. Correspondingly, 17(c) = it(eic) = ^(c 1 ) where u 
is a utility function of a scalar argument. Our presentation is oriented to the 
power utility function u 7 (x) = x 7 /y with 7 g] 0,1[. The case of 7 < 0, where, by 
convention, uo(x) = lnx, is of interest but it is not covered by the present study. 

Remark 2. We consider here a model with mixed “regular-singular” controls. In 
fact, the assumption that the consumption process has an intensity c = (ct) and the 
agent’s utility depends on this intensity is not very satisfactory from the economical 
point of view. One can consider models with an intertemporal substitution and 
the consumption by “gulps”, i.e. dealing with “singular” controls of the class A x 
and the goal functionals like 


Jt := / e~ 0s U(C s )ds , 


where 


C s = [ K(s, r)dC r 

Jo 


with a suitable kernel K(s,r ) (the exponential kernel e 7 ^ s is the common 
choice). 


4 Continuity of the Bellman Function 

Proposition 4.1 Suppose that W(x) < 00 for all x G int K. Then W is continuous 
on int K . 

Proof. First, we show that the function W is upper semicontinuous on int A'. Sup¬ 
pose that this is not the case and there is a sequence x n converging to some 
xo G int A' such that limsup ra W(x n ) > W(xo). Without loss of generality we way 
assume that the sequence W{x n ) converges. The points % = (1 + l/k)xo, k > 1, 
belong to the ray R+*o and converges to xo■ We find a subsequence x nk such that 
Xk —K Xn k for all k > 1. Indeed, since 

X}. = (1 + l/fc)* 0 G xq + int A', 

there exists e*. > 0 such that 

Xk + Oe k (0) G *0 + int K. 


It follows that 

Xk + (Xn - xo) + Oe k (0) G Xn + int K 

and, therefore, Xk G x n + int A' for all n such that \x n — *o| < £fc- Any strictly 
increasing sequence of indices nk with \x nk — xo\ < gives us in a subsequence of 
points x nk having the needed property. The function W is increasing with respect 
to the partial ordering >k- Thus, 

lim W{xk) > \ivoW{xn k ) > W(xo). 
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On the other hand, the function A h4 W(\xq) is right-continuous at A = 1 and, 
hence, lim*. W(xk) = W(xo). This contradiction shows that W is upper semicon- 
tinuous on int K. 

Let us show now that liminf n W(x n ) > W(x o) as x n —> xo, i.e. W is lower 
semicontinuous on int K. 

Fix e > 0. Due to the finiteness of the Bellman function there are a strategy n 
and T £ R+ such that for 9 = 9 X ° ,7r we have the bound 

pT A 6 

E / e~ fis U(c s )ds > W{x o) - e. 

Jo 

It remains to show that 


liminf 9 n A T > 6 A T a.s., (12) 

n 

where we use the abbreviation 9 n '■ = 9 Xn ’ n . Indeed, with this bound we get, using 
the Fatou lemma, that 

p6 n AT p6 n AT 

liminf W(x n ) > liminf E / e~^ s U(c s )ds > E liminf / e~^ s U(c s )ds 
n n J Q n Jo 

pOAT 

>E e~ Ps U{c s )ds > W(xo) - e 

Jo 

and the claim follows since e is arbitrarily small. 

To prove (1121) . we observe that by virtue of m on the interval [0, 9 n A 9 A T] 
we have the representation 

V) Xn ’ 7r - V t X0 ’ n = D Xn - X0 S t 

implying that 

sup \V t x ^ - V t Xo,7T \ < \x n - *o|, 

t< 6 >„A 6 lAT 

where S := sup t<T l^tl- Fix arbitrary, “small”, <5 > 0. For almost all uj the distance 
p(ui) of the trajectory V Xo,7r (uj) from the boundary dK on the interval [0,0(w) A 
T — S\ is strictly positive. The above bound shows that for sufficiently large n the 
trajectory y Xn, 7 r (w) does not deviate from more than on p(u >)/2 on the 

interval [0, 9 n (uj) A 6(u>) AT]. It follows that 9 n (oj) > 9(uj) A T — S. Thus, 

lim inf 9 n AT>9AT — S a.s. 

n 

and (m holds. □ 


5 The Hamilton—Jacobi—Bellman Equation 

Let G := (— K ) n dO\(0) where O r {y ) := {x £ R d : \x — y\ < r}. The set G is 
compact and —K = con eG. We denote by Eq the support function of G, given by 
the relation Eq{p) = sup xgG pa;. The convex function U*(.) is the Fenchel dual of 
the convex function —[/(—.) whose domain is —C, i.e. 

U*(p) = sup(17(a;) — px). 

X^lC 
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We denote by Ci(A') the subspace of the space of continuous functions / on 
A such that sup^g^g \f(x)\/(l + |x|) < oo. In other words, Ci(A) is the space of 
continuous functions on K of sublinear growth. The notation / £ C 2 (x ) means 
that / is smooth in some neighborhood of x. 

Let / £ Ci (A) n C 2 (int A'). Using the abbreviation 

I(z,x) ■ I{z: X+D x z£int K} -^int k( x T F) x Z ) 
we introduce the function 

l(f,x) := j [( f(x + D x z)I(z,x)-f(x))-D x zf'(x)]n(dz ), a; £ int A'. 

It is well-defined and continuous in x. Indeed, fix xo £ int K. Let e £]0,1] be such 
that the ball O 4 £ (xo) C K. With this choice x + D x z £ C> 2 e(*o) when x £ O e (x o) 
and \z\ < 5 := e/(l + |xo|)- Using the Taylor formula for such value of a and the 
sublinear growth of / for 2 with \z\ > S we obtain the following uniform bound for 
x £ O e {x 0 ): 

\(f(x + D x z)I(z,x) - f(x) - Dx.zf'(x )| < m\z\ 2 1 0s{xo) (z) + K 2 \z\I O c {0 )(z). 

It implies the needed integrability and the continuity of the integral in x. 

We introduce a function of five variables by putting 

F(X,p,X(f, x),W, x):= max{ F 0 (X,p,X(f , as) ,W,x) + U* ( p ), E G (p)}, 

where X belongs to Sd, the set of d x d symmetric matrices, p, x £ R c/ , W £ R, 
/ £ Ci (A') n C 2 (x), and the function Ao is given by 

F 0 (X,p,X(f,x),W,x) := ±tiA(x)X + p(x)p + X(f,x)-0W 

where A{x) is the matrix with A lJ (x) := a ZJ x l x J , and p(x) is the vector with 
components jA{x) := p l x l , 1 < i,j < d. 

In a more detailed form we have that 

d d 

F 0 (X,p,X(f,x),W,x) = - ^2 a %0 x l x J X %3 + ^ iix % p l + X (/, x) - ftW. 

i,j =1 i= 1 

Note that Fo is increasing in the argument / in the same sense as X. 

If 0 is a smooth function, we put 

C<j)(x) := F(<j>"(x),(l>'(x),T((l>,x),<l>{x),x). 

In a similar way, £0 corresponds to the function Fq. 

We show, under mild hypotheses, that IT is a viscosity solution of the Dirichlet 
problem for the HJB equation 

F(W"(x),W'(x),X(W,x),W(x),x) = 0, a; £ int A, (13) 

W{x) = 0, x £ dK, (14) 

with the boundary condition understood in the usual classical sense and establish 
a uniqueness result for this problem. 
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6 Viscosity Solutions for Integro-Differential Operators 

Since, in general, W may have no derivatives at some points x £ intA - (and this 
is, indeed, the case for the model considered here), the notation (1131) needs to be 
interpreted. The idea of viscosity solutions is to substitute W in F by suitable 
test functions. Formal definitions (adapted to the case we are interested in) are as 
follows. 

A function v £ C'(A') is called viscosity supersolution of (1131) if for every x £ int A 
and every / £ Ci(A') D C 2 (x) such that v(x) = f{x) and v > f the inequality 
Cf{x) < 0 holds. 

A function v £ C(K) is called viscosity subsolution of (1131) if for every x £ int A 
and every / £ C'i(A) D C 2 (x) such that v{x) = f(x ) and v < f the inequality 
Cf{x) > 0 holds. 

A function v £ C(K) is a viscosity solution of (1131) if v is simultaneously a 
viscosity super- and subsolution. 

At last, a function v £ C\{K) DC 2 (int A) is called classical supersolution of (1131) 
if Cv < 0 on int A. We add the adjective strict when Cv < 0 on the set int A. 

For the sake of simplicity and having in mind the specific case we shall work 
on, we incorporated in the definitions the requirement that the viscosity super- 
and subsolutions are continuous on K including the boundary. For other cases this 
might be too restrictive and more general and flexible formulations can be used. 

Lemma 6.1 Suppose that the function v is a viscosity solution of m- If v is twice 
differentiable at xq £ int A, then it satisfies m at this point in the classical sense. 

Proof. One needs to be more precise with definitions since it is not assumed that 
v' is defined at every point of a neighborhood of xq. “Twice differentiable” means 
here that the Taylor formula at xo holds: 


v(x) = P 2 (x - xq) + (x - xo) 2 h(\x - a:o|) 


where 



r > 0. Note that O r j\ Xa \ C A r . Hence, I7(rf) < oo. 

Let e £]0,1]. We choose a number <5o €]0,1[ such that so + £Ji o (0) ^ int A and 
|h(s)| < £ for s < So- Put 8 := <5o/(l + |*o|)- Take A £]<5,<So[ sufficiently close to 8 
to insure that xq + C>zl(0) C int K and n{P/^ \ Pf) < e. 

We define the function f e £ C'i(A') D C 2 (x o) by the formula 



P 2 (x - xo ) + s(x - xo) 2 , 

g(x) V v(x), 

v(x), 


x £ x 0 + 0,5(0), 
x € xo + Oa(0) \ 0^(0) 
X £ Xq T 0^(0), 



where 
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Clearly, f s (x o) = v(xo) and f s > v. Since v is a viscosity subsolution, £f s (x o) > 0. 
Note that 

n 

\£fe(xo) - Cv(x o)| < e'^a\x o) 2 + T(f e -v,x 0 ), 

i =1 

with 

T(fe-v,x o) = J(f E - v){x 0 + D Xo z)I {xo+DxoZeintK} n(dz). 

Let us check that I(/ e — v,xo) < ks. Indeed, 

(fs ~ v )( x o + Ae 0 2 ) < e (Dx 0 z) 2 Ir 5 + MIr^\r s 

< emm{\x 0 \ 2 \z\ 2 ,5 2 } + MIr A \r s - 

where A/ = 1 + sup, yeC ) i( - 0 ) l-Fb(j/) |- It follows that 

T((/e - v)I 0 ^(o),xo) < e(l + \x 0 \) 2 I k| 2 A 1 II(dz) + Me. 

Letting e tend to zero, we obtain that Cv(x 0 ) > 0. Arguing similarly with e < 0, 
we get the opposite inequality. □ 


7 Jets 

Let / and g be functions defined in a neighborhood of zero. We shall write /(.) 

g(-) if f(h) < g(h ) + o(\h\ 2 ) as \h\ -> 0. The notations /(.) > g(.) and /(.) « p(.) 

have the obvious meaning. 

For p € R d and A € 5^ we consider the quadratic function 

Qp,x(z) :=pz + (l/2){Xz,z), z € R d , 

and define the super- and subjets of a function v at the point x\ 

J + v(x) := {(p, X) : v(x + .) < v(x) + Q p ,x(-)}, 

j~v(x) := {(p,x) : v(x + .) > v(x) + Q Pl x(-)}- 

In other words, J + v(x) (resp. J~v(x)) is the family of coefficients of quadratic 
functions v(x) + Q p ,x(y— •) dominating the function v(.) (resp., dominated by this 
function) in a neighborhood of the point x with precision up to the second order 
included and coinciding with v(.) at this point. 

In the classical theory developed for differential equations the notion of vis¬ 
cosity solutions admits an equivalent formulation in terms of super- and subjets. 
Since the latter are “local” concepts, such a characterization is not possible for 
integro-differential operators. Nevertheless, one can construct from semijets test 
functions with useful properties. 

The following lemma claims for v G C\(K) with any element (p, X) G J + v(x) 
x G int K one can relate a test function dominating v, arbitrary close to v in the 
uniform metric, touching v at the point x, smooth at a neighborhood of x and 
having at this point the first and the second derivatives coinciding with p and A'. 
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Lemma 7.1 Let v £ Ci(K) and let a > 0. Let x £ int K and let ( p,X ) £ J + v(x). 
Then there exist a number ao £]0, 1[ and a C 2 -function r : R d —t R with compact 
support such that 

lim \u\~ 2 r(h) = 0, (15) 

|/t|—s-o 

and the function fo : K —» R given by the formula 

fo(x + h):= [(v(a) + Q p ,x{h) + r(/i)) Vt(x + /i)] A(v(x + h) + a), x + heK, (16) 

has the following properties: 

fo(x + h) = v(x) + Q p ,x(h) + r(h), h € O ao { 0), 

v < fo < v + a on K, f 0 (x) = v(x), f' 0 (x) = p, fo (x) = X. 

In particular, if v is a subsolution of the HJB equation, then Cf < 0 on int A'. 

Proof. Take ao e]0,1[ such that the ball O 2 a 0 {x) = {y € R d : |y — x\ < 2 a} lays in 
the interior of K. By definition, 

v{x + h) - v(x) - Q p ,x(h) < \h\ 2 g>(\h\), 

where <p(u) —> 0 as it 4, 0. We consider on ]0,ao[ the function 

S(u) := sup —! ,2 («(* + h) - v(x) - Q p ,x(h)) + < sup V + (y). 

{h: |h|<u} l^l { y : 0<y<w} 

Obviously, S is continuous, increasing and 5(u) —> 0 as u f 0. We extend 5 to a 
continuous function on R+ with 8(u) = 0 for u > 1. 

The function 

q i-2u i>2t] 

^(m) : = 3 f J <5(C)dC^ 

vanishes at zero with its two right derivatives; u 2 5{u ) < A(u) < M 2 5(4tt). It follows 
that the function r : h i-i- A(\h\) has a compact support, belongs to C 2 (O ao {0 )), 
its Hessian vanishes at zero, and 

v(x + h) - v(x) - Q p ,x{h) < \h\ 2 S(\h\) < A(\h\) = r(h), h £ O ao ( 0) 

Thus, the function y v{x) + Q Pt x{y — x) + r(y — x) dominates v on the ball 
O ao (x). Without loss of generality, diminishing ao if necessary, we may assume 
that it is dominated by v + a on this ball. Now the assertion of the lemma is 
obvious. □ 

The corresponding assertion for J~v(x) also holds — with obvious changes in 
the formulation. 

For the proof of the uniqueness theorem we need specific families of test func¬ 
tions coinciding with sub- and supersolutions outside a neighborhood of x. To this 
end we introduce the following definitions. 

Let 0 < a < a! . We say that a continuous mapping £ a a / : R d —> [0,1] is an 
{a,a')-cutoff function if £ a>a / = 1 on O a ( 0) and £ a a i = 0 outside C> a /(0). If L is a 
linear subspace of R d we define the cylindrical (a, a')-cutoff function a , by putting 
£a,a'{ x ) = (,a,a'{PLx) where Pl is a projection of x onto L. 
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It is clear that in the notation of the above lemma for any a' €]0, ao] the 
functions f : K —¥ R given by the formulae 

f(x + h) := (v{x) + Q P 'X(h) + r(h))£ a!a >(h) + v(x + h)(l - £ a ,a'(h)), x + h 
f(x + h) := fo(x + h)£a, a '(h) + v(x + h)(l -£a >a '(h)), a c + heK, (18) 

will satisfy all the properties claimed for /o. 

The following lemma will be used in the case when D = D x = diaga;, D = 
D y = diagy, and x, y has no zero components. 

Lemma 7.2 Let D, D be two invertible linear operators on R d . Let £ 0 a / be a ( a,a’)~ 
cutoff function. Then there is a {a, a!)-cutoff function £a,a' with arbitrary small a < 
a.|1 1 1| 1 and arbitrary a' > a'\\DD~ 1 \\ such that 

la,a>(Dz) = £ a y(Dz) V Z £ R d . 

Proof. Put 

la.a'(«) := « C R d . 

Then Ca,a'(w) = 1 if |u| < a||I)H _1 |p 1 , and £a,a'{ u ) = 0 if > a'. The last 

inequality holds when |m| > a , ||(.D.D -1 )~ 1 || = a'\\DD~ 1 \\. □ 

Remark 7.1 The assertion of the lemma remains true for cylindrical cut-off func¬ 
tions in the situation where D and D are two symmetric operators with the com¬ 
mon image space L = ImB = ImD. The norm in formulation is the norm in L of 
their restrictions. 


8 Supersolutions and Properties of the Bellman Function 

8.1 When is the Bellman Function W Finite on K7 

First, we present sufficient conditions ensuring that the Bellman function W of 
the considered maximization problem is finite. 

Functions we are interested in are defined in the solvency cone K while the 
process V may jump out of the latter. In order to be able to apply later the Ito 
formula we stop V = V x '^ at the moment immediately preceding the ruin and 
define the process 

v = v 8 - = VI [0 , e[ + v s _i [6t00[ , 

where 9 is the exit time of V from the interior of the solvency cone K. This process 
coincides with V on [0. 6[ but, in contrast to the latter, either always remains in K 
(due to the stopping at 9 if V@ — € int K) or exits to the boundary in a continuous 
way and stays on it at the exit point. 

Since Vt = V tA g — we obtain from the equation (J8]) the represen¬ 

tation 

/■tAf? rt/\6 n 

Vt=x+ / Dy (fids + Sdw s ) + / / Dy z(p(dz,ds) - q(dz,ds)) 

Jo s Jo J B 

-AV 0 I [ejOo[ (t) + B t -C t . 
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Let <P be the set of continuous functions f : K —¥ R+ increasing with respect to 
the partial ordering >x and such that for every x £ int K and 7 r £ AS the positive 
process X f = X given by the formula 

X t f := e-*V(Vi)/[o,9[(t) + Jt (19) 

is a supermartingale. 

The set $ of / with this property is convex and stable under the operation A 
(recall that the minimum of two supermartingales is a supermartingale). Any con¬ 
tinuous function which is a monotone limit (increasing or decreasing) of functions 
from also belongs to <P. 

The interest to the processes X? with / £ is explained by the following: 
Lemma 8.1 (a) If f £ <&, then W < f. 

( 6 ) Let y £ dK. Suppose that for every e > 0 there exists f e £ such that f e ( y ) < e. 
Then W is continuous at y and W(y) = 0. 

Proof, (a) On the boundary dK the inequality is trivial. Using the positivity of /, 
the supermartingale property of X\ and, finally, the monotonicity of / we get, 
for £ int K, the following chain of inequalities leading to the required property: 

EJt < EX{ < /(Vo) = /(Vo) < /(Vo-) = /(*)• 

(6) The continuity of the function W at the point y £ dK follows from the 
inequalities 0 < W < fe- □ 

Remark. Recall that Proposition 14.11 assets that the function W, if finite, is con¬ 
tinuous on the interior of K. Thus, the above lemma implies that W is continuous 
on int K if <P is not empty. If 4> is reach enough to apply ( b ) at every point of the 
boundary, then W is continuous on K and vanishes on the boundary. 

Lemma 8.2 Let f : K — » R+ be a function in C\ (K) fl C 2 (int A'). Iff is a classical 
supersolution of m, then f £ L>, i.e. f is increasing with respect to the partial ordering 
>K and X^ is a supermartingale. 

Proof. First, notice that a classical supersolution is increasing with respect to the 
partial ordering >k- Indeed, by the finite increments formula we have that for any 
x, h £ int K 

f(x + h) - f(x) = f'{x + dh)h 

for some •& £ [0,1]. The right-hand side is greater or equal to zero because for 
the supersolution / we have the inequality SQ(f'(y)) < 0 whatever is y £ int A, 
or, equivalently, f{y)h > 0 for every h £ K, just by the definition of the support 
function Eq and the choice of G as a generator of the cone —K. By continuity, 
f{x + h) — f(x) > 0 for every x, h £ K. 

Let 6 n := inf{t: dist(Vt, OK) < 1/n}. The stopped processes V evolves in 
int K. Thus, we can apply the “standard” Ito formula to e~^ t f(Vt) and obtain, for 
t <6, that 

e~ 0t f(V t ) = f(x) + [ e- ps f\V s -)dVs-p [ e- ps f(V s -)ds 
Jo Jo 

+ \ I e~P s trA(V 3 —)f"(V a —)ds 
+ ^e _/ 3 s [/(V s - +AV S )~ f(V s -) - 

S<t 


}'{V s -)AVs). 
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Taking into account that the processes Y and B do not jump simultaneously 
and the ruin cannot happen due to a jump of B we get that 


J2 e ~ Ps lftV- + A V t ) - /(V._) - f'(V s -)AV s ] - e-P e f'(Ve_)AV e I {g} (t) 

S<t 

-e-P e f(Ve-)I {e} (t) 

= + AV s )I intK (y s - + AV a ) - f(\ 4-) - f'(V a -)AV s ] 

S<t 

= f f e~ Ps [f(y s - + D Vs _z)I(V s -,z) - f(V a -) - f'(V 3 -)D Vs _z)]I {ABs=0} p(ds,dz) 
+ X> _/3s [/04- +AB S ) -f(V s -) - f{V s -)AB s ] 

S<t 


I J e l3S ['"] I { AB <‘=°}( p ( ds ’ dz ') ~ n ( dz ) ds ) + J f < 


+ + AB.) - f(V B -) - f(Vs-)AB a ] 

S<t 


[...]II(dz)ds 


where we replace in the integrals by dots the lengthy expression 
f(Ve- + Dy s _z)I(V s -,z) - f{Vs-) - f{Vs-)Dy s _z. 


Noting that 

X{ = e~ 0t f(V t ) - e~^f(V e _)I {e} {t) + J T 

and using the equation m and the above formulae we obtain, after regrouping 
terms, the following representation 


X{ = 


rtAO 

/(*)+ / 

Jo 


e-^lCofiVs) - c s f\Vs) + U{c s )]ds + R t + m t , (20) 


where 


Rt ■■= 


rtA0 

■= / e~ Pa f{V s -)dB c s + Ve" /3s [/(t>_ + AB S ) - f{V s -)] (21) 

% '° S<t 


and m is the local martingale 

rtAO 


m t = 


/ 

Jo 


s f(Vs-)Dy ~dw s 


rtAd r 

Jo ./' 


3 [/(Vs_ + Dy 3 _ z)I(Vs -, z) - /(Vs-)] (p(dz, ds) - n(dz)ds). (22) 
By definition of a supersolution, for any x £ int K, 

C 0 f(x) < -U*(f'(x)) < cf\x) - U(c) Vc € C. 


Thus, the integral in (1201) is a decreasing process. The process R is also decreasing. 
Indeed, the terms of the sum in (ED are less or equal to zero in virtue of the 
monotonicity of / and 

f(V s -)dB c s = I {ABs=0} f'(V s -)B s d\\B\\ s 
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where f'(V s -)B s < 0 since B takes values in — K. Let a n be a localizing sequence 
for m. Taking into account that Xf > 0, we obtain from m that for each n 
the negative decreasing process RtAa n dominates an integrable process and so it 
is integrable. The same conclusion holds for the stopped integral. Being a sum of 
an integrable decreasing process and a martingale, the process x( Acrn is a positive 
supermartingale and, hence, by the Fatou lemma, X? is a supermartingale as well. 
□ 

Lemma [8.21 implies that the existence of a smooth positive supersolution / of 
m ensures the finiteness of W on K. We discuss a method how to construct 
supersolutions in Section 12. 

Remark. Let O be the closure of an open subset O of K and let / : O — > R+ 
be a classical supersolution in O increasing with respect to the partial ordering 
>K- Let x G O and let r be the exit time of the process V X '~ K from O. The above 
arguments imply that the process x( A is a supermartingale and, therefore, 

E[e-^ tAT) f(VtAT)I l0 , e[ (t A r) + JtAr] < fix). (23) 


8.2 Strict Local Supersolutions 

For the strict supersolution we can get a more precise result which will play the 
crucial role in deducing from the Dynamic Programming Principle the property 
of IT to be a subsolution of the HJB equation. 

Fix x G int I\ and a ball O r (x) C int K such that the larger ball & 2 r[x) C int K. 
We define t 77 = r * as the exit time of V^’ x from O r (x), i.e. 

t 77 := inf{t > 0 : \V^’ X — x\ > r}. 

Lemma 8.3 Let f G Ci(K) DC 2 (C> 2 r(x)) be such that Cf < —e < 0 on O r {x). Then 
there exist a constant = r/ e > 0 and an interval ] 0 , to] such that 

sup EX{^ < f{x) -r]t V t G]0, t 0 |u 

Proof. We hx a strategy n and omit its symbol in the notations below. In what 
follows, only the behavior of the processes on [0,r] does matter. Note that \V T — x\ > 
r on the set (r < oo} and r < 6. As in the proof of Lemma 18.21 we apply the Ito 
formula and obtain, with the same notations, m and (1221) . the representation 

X t f AT := e~^ tAT) f(V t Ar)I [0A (t A r) + J& T 

rt/\T 

= f{x)+ / e~ 0s {C o f + U*)(V s )ds 
Jo 

rtAr 

- / e- 0s [U*{Vs) + CsfiVs) - U(c s )}ds + RtAr + mtAT- 

■to 

Due to the monotonicity of / we may assume without loss of generality that on 
the interval [0, r] the increment ABt does not exceed the distance of V s - to the 
boundary of O r (x). In other words, if the exit from the ball is due to an action 
(and not because of a jump of the price process), we can replace this action by 
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a less expensive one, with the jump of the process V in the same direction but a 
smaller one, ending on the boundary of the ball. So, \ABt\ < 2 r for t < r. 

By assumption, for y £ O r (x ) we have the bounds (£o f+U*)(y) < —s (implying 
that the first integral in the right-hand side above is dominated by — e (tAr)) and 
E'df'iy )) < —e. The latter inequality means that kf' {y) < — e\k\ for every k £ —K 
(therefore, we have the inclusion f'(O r (x)) C int K*). In particular, for s £ [0,r] 

f{V s -)B s < -e\B a \, [/( Vs- + AB S ) - f(V B -)] < -e\AB s \. 

Since \V S - — x\ < r for s £ [0,r], we obtain, using the finite increment formula and 
the linear growth of /, the bounds 

\f(V s - + D Vs _z) - /(V s _)] 2 7(V s _^)7 { |,|< 1/2} < k|^| 2 7 { | 2 |< 1/2} , 

If (V s - + Dv s _z) - f(V s -)\I(V s -,z)I {lz]>1/2} < k{ 1 + \z\)I m>1/2} , 
and, as I(V s -,z) = 1 when \z\ < r/(\x\ +r), 

f(V s -)(l-I(V s -,z))<Kl { \ z \> r/2} 

for some constant k independent on the strategy. Thus, the integrand in the 
stochastic integral with respect to the centered Poisson measure in (1221) for t < r 
is bounded by the function \z\ 2 A \z\ multiplied by a constant while the integrand 
in the integral with respect to the Wiener process is bounded. It follows that the 
local martingale (mtAr)t>o is a martingale and Emt/\ T = 0 . 

The above observations imply the inequality 

EX{& < f{x) - e~ 0t EN t , 

where 

/•tAr ntAr 

N t :=£(t/\r)+ / H(c s ,f\V s ))ds + e |5 s |d||B|| s 
Jo Jo 

with 77(c,p) := U*{p) +pc— U{c) > 0. It remains to verify that ENt dominates, on 
a certain interval ] 0 ,to], a strictly increasing linear function which is independent 

of 7T. 

The process Nt looks a bit complicated but we can replace it by another one 
of a simpler structure. To this end, note that there is a constant k (“large”, for 
convenience, k > 1 ) such that 

inf 77(c,p) > ^|c|, Vc£C, |c| > k. 

pe/'(Or( x)) 2 

Indeed, being the image of a closed ball under continuous mapping, the set f'(O r (x)) 
is a compact in int K*. The lower bound of the continuous function U* on f'(O r (x)) 
is finite. For any p from f'(O r (x)) and c £ C C K we have the inequality pc/|c| > e. 
At last, U{c)/\c\ —> 0 as c —» oo. Combining these facts we infer the claimed in¬ 
equality. Thus, for the first integral in the definition of Nt we have the bound 

/•tAr ptAr 

J ^ H(c s j'(Vs))ds>- I 


7 {|Ca|>«}l Cs l dS - 
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The second integral in the definition dominates Ki||-B||tAr for some ki > 0. To 
see this, let us consider the absolute norm |.|i in R d . In contrast with the total 
variation ||_B|| which is calculated with respect to the Euclidean norm |.|, the total 
variation of B with respect to the absolute norm admits a simpler expression 
Var B l where Var B l is the total variation of the scalar process B l . Obviously, 


i*ix = 5>‘i = E 

i i 


dB l 


dB l 

dVar B i 

d E l VarS l 

d\\B\\ 

i 

dVar B 1 

d\\B\\ 

d\\B\\ 


But all norms in R d are equivalent, i.e. k _1 |.| < |.|i < «|-| for some strictly pos¬ 
itive constant k. The same inequalities relate the corresponding total variation 
processes. The claimed property follows from here with the constant k\ = k -2 . 

Summarizing, we conclude that it is sufficient to check the domination property 
for ENt with 

rtAr 

N t :=tAr+ / I{\ Cs \> K }\ c s\ds+ \\B\\ tAr . (24) 

Jo 

These processes N = N n have a transparent dependence on the control. The idea 
of the concluding reasoning is very simple: on a certain set of strictly positive 
probability, where one may neglect the random fluctuations, either r is “large”, 
or the total variation of the control is “large”: one can accelerate exit only by an 
intensive trading or consumption. 

The formal arguments are as follows. Using the stochastic Cauchy formula 0 
and the fact that So +(Y l ) = So(Y l ) = 1, we get immediately that there exist a 
number to > 0 and a measurable set B with P(B) > 0 on which 

\V x ’ n - *1 < r /2 + 2(||B|| + He'll) on [0,i 0 ] 


whatever is the control 7 r = ( B,C ). Of course, diminishing to, we may assume 
without loss of generality that Kto < r/8. For any t < to we have on the set 
r Pi {t < t} the inequality ||B|| T + ||C|| T > r/4 and, hence, 

Nt > ||R||r + ||C||r - I -f{|c s |<K}l Cs l ds > J - Kt 0 > Kto >to>t. 

On the set B C\ {t > t} the inequality N't > t is obvious. Thus, ENt > tP(T) on 
[ 0 , to] and the result is proven. □ 


9 Dynamic Programming Principle 

The aim of this section is to establish the following two assertions which will serve 
to derive the HJB equation for the Bellman function. For the considered model, 
they constitute an analog of the classical Dynamic Programming Principle. The 
latter is usually written in the form of a single identity (see the remark at the end 
of the section), but for our purpose this form, more precise, is needed. 

Lemma 9.1 Let Tf be the sets of finite stopping times. Then 

W(x) < sup inf E(j^ + e~ 0 T W(Vr' 7 T )I, T<g} ) . (25) 

Tre^s re7> v *) 
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Lemma 9.2 Suppose that W is continuous on int K. Then for any r £ 7/ 

W(x) > sup£(j T T + e~ 0T W{V x ^ )I {T< o}) ■ (26) 

We work on the canonical filtered space of cadlag functions equipped with the 
measure P which is the distribution of the driving Levy process. The generic point 
uo = cl), of this space is a d-dimensional cadlag function on R + , zero at the origin. 
Let Pf := <t{oj s , s < t} and Pt '■= n £ >od r t ° +£ . We add the superscript P to denote 
(j-algebras augmented by all P-null sets from 17. Recall that P t ' coincides with 
pf (this assertion follows easily from the predictable representation theorem). The 
Skorohod metric makes 17 a Polish space and its Borel a-algebra coincides with 
Poo, for details see [?]. 

Since elements of 17 are paths, we can define such operators as the stopping 
ui. H > ojf, s > 0, where uif = io s /\. and the translation u>. t->- u s -|-. — oj s . Taking Doob’s 
theorem into account, one can describe J^-measurable random variables as those 
of the form g(w.) = g(w s ) where g is a measurable function on 17. 

We define also the “concatenation” operator as the measurable mapping 


g : R + x 17 x 17 —> 17 


with gt(s,u].,u>.) — (jJtI[o,s[{t) + (pt—s + w s)7[ Sj0 o[(f)- 

Notice that 

gt(s,U0 S ,U0. + s - LOs) = UJ t - 
Thus, 7r(cu) = Tv(g(s,oj s ,oj. +s -ui s )). 

Let 7r be a fixed strategy from and let 9 = 9 X,7T be the exit time from int K 
for the process V x,7r . 

Recall the following general fact on regular conditional distributions. 

Let £ and g be two random variables taking values in Polish spaces X and Y 
ec l u ipped with their Borel a-algebras X and y. Then £ admits a regular conditional 
distribution given p = y which we shall denote by p^ v (r,y). This means that 
p^l^(.,y) is a probability measure on X, p^^P ,.) is a T-measurable function, and 


E(m,v) 1 * 7 ) 


I f(x,y)Pt;\ v (dx,y) 


y=T] 


{a.s.) 


for any X x T-measurable function f{x,y ) > 0. 

We shall apply the above relation to the random variables £ = ( uj.+ t — u> T ) 
and y = ( t,u > t ). It is well-known that the Levy process starts afresh at stopping 
times, i.e. the measure P(.) itself (not depending on y) is the regular conditional 
distribution p^ v {.,y). 

At last, for fixed s and w s , the shifted control 7 r._|_ s (gi(s,u; s ,a).)) is admissible 
for the initial condition V x,7r (uj) when s < 9{uj). Here we denote by u>. a generic 
point of the canonical space. 

Proof of Lemma \M For arbitrary n G Aa and 7/ we have that 

POO 

EJZo = EJr + Ee~ PT I {T<e} / e~^ r U(c r + T )dr 

Jo 

= EJ? + Ee-P T I {T<e} E( I e _/ 5 r l/(c r + r )dr|(T ) w T )). 




20 


Dimitri De Valliere et al. 


According to the above discussion we can rewrite the second term of the right-hand 
side as 

Ee~ l3r I{ T<e } j (^j e^' 3r t/(c r+T (p(r,a; r ,w)))drj P(dw) 

and dominate it by Ee~^ T I^ T<g ^W(Vr’ n )- Thus, 

EJZo < EJ? + Ee- pT I {T<e} W(V^). 

This bound leads directly to the announced inequality. □ 

Proof of Lemma El Fix e > 0. By hypothesis, the function W is continuous on 
int A. For each x G int A we can find an open ball O r (x ) = x + O r { 0) with 
r = r(e,x) < £ contained in the open set {y G intA : \W(y) — W(x)\ < e}. 
Moreover, we can find a smaller ball Of(x) contained in the set y(x) + K with 
some y(x) G O r {x). Indeed, take an arbitrary xq G int A'. Then, for some 5 > 0, 
the ball xo + 0,5(0) C K. Since A is a cone, Xxo + 0*5(0) C A for every A > 0 and 
this inclusion implies that 


x + 0*5(0) C x - Xx 0 + K 

Clearly, the requirement is met for y(x) = x — Xxo and f = XS when A|xo| < r and 
A 5 < r. The family of sets O f ( x y 2 (x), x G int A, is an open covering of int A. But 
any open covering of a separable metric space contains a countable subcovering 
(this is the Lindelof property; in our case, where int K is a countable union of 
compacts, it is obvious). Take a countable subcovering indexed by points x n - For 
notational simplicity, we shall denote the open balls 0f(x„)/ 2( x n) by On and y(x n ) 
by y n . 

Let 7r n = (B n , C n ) G Aa n be an e-optimal strategy for the initial point y n , i.e. 
such that 

Aj£r > w(y„) - £- 

Let 7r G Aa be an arbitrary strategy. Put 

p := inf{j > 1 : V T X ’ W € Ojj, 

Let us introduce the strategy 

^ ■ ^f[0,r] "F (0;0)f]r,oo[ 

and the predictable stopping times 77 . := t + 1/k. Finally, put 

OO 

F ■ — kI[ 0 ,t] + \{yn — Vr k ,0) + 7r ]I[Tk, (> o[I{p= n }I{y*^' _y n £K}I{ T k < 8} 

n= 1 

where n n ’ k is the translation of the strategy n n : namely, for a point u>. with t(uj.) = 
s < 00 we have 

At ’ t (w.) := 7T t.-s-l/ki^.+s+l/k — w s+l/fe)' t — Sl/k. 

In other words, the strategy if coincides with 7r on [0,r[, is zero on the interval 
[r, Tfc[ and with the shift of -K n on [77., 00[ when V)? ,7r is in O n and V ^. ,7r — y n G A; 
the correction term guarantees that in the latter case the trajectory of the control 
system corresponding to the control n passes at time 77, through the point y n . 
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One can check that nI^ o gx ,^ £ A x . 

Now, using the same considerations as in the previous lemma, we have: 

00 POO 

W(x) > EJI =EJ? + Y j Vhp^hrKe}!^'_ yneK} / e~P s U{c n s )ds 

n= 1 k *' T fc 

00 

>EJ? + J2 EI {p=n} I {T<e} e- pT (W(y n ) - e) 

71=1 

> EJ? + Ee~ pT W{V^)I {T<e} - 2e. 

Since 7r and e are arbitrary, the result follows. □ 

Remark. The previous lemmas imply that for any r £ 7/ the following identity 
holds: 

W(aO = sup E (j? +e-P T W(V?’ n )I {T<e} ) . 

It can be considered as a form of the dynamic programming principle but, seem¬ 
ingly, it is not sufficient for our derivation of the HJB equation. 


10 The Bellman Function and the HJB Equation 

Theorem 10.1 Assume that the Bellman function W is in C{K). Then W is a vis¬ 
cosity solution of El- 

Proof. The claim follows from the two lemmas below. □ 

Lemma 10.2 If W is in C(int A') then W > 0 is a viscosity supersolution of 11 S\i . 

Proof. Let x £ int I\ and let £ C' 1 (A')nC 2 (a;) be a function such that 4>(x) = W(x) 
and W > 4> on K. 

Fix an arbitrary point m £ K. Let e > 0 be sufficiently small to guarantee that 
x — em £ O r (x). The function W is increasing with respect to the partial ordering 
generated by K. Thus, 

4>{x) = W{x) > W(x — em) > 4>{x — em). 

Taking a limit as e —t 0 in the inequality e~ 1 {(j>{x — em) — 4>{ x )) — 0, we obtain 
that —m4>'(x) < 0 and, hence, Eq(4>'(x)) < 0. 

Take now n with Bt = 0 and ct = c £ C for all t. Let r r = t? < 9 be the 
exit time of the process V = V x,w from the ball O r (x)\ obviously, r r < 6. The 
properties of the test function and the inequality ([261) imply that 

H x ) = w(x) > E (j? ATr + e -/3(*A ^W(VtAT r )l { tAT r< e}) 

> E (jTat, + e-^ tATr U(Vt ATr )I{tAT r <e}) ■ 
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We get from here, using the Ito formula (1201) . that 


0 > E 


(/ e l3s U(c s )ds + e p{tATr) <f>( v t/\T r )I{tAT r <e}^ ~ H x ) 

ptAr r 

Tr<e} / e- ps [Co<t>{Vs)-c^(y a ) + U{c)]ds 

JO 

> min [C 0 <f>(y)-af>'(y) + U(c)]EI[ t ^ T <g - i ^ (l - e ^ /j(tATr ' > ) 
yeOr(x) P ' ' 


EI{tA 


Dividing the resulting inequality by t and taking successively the limits as t and r 
converge to zero we infer that Co<j>{x) — aj>’(x) + 17(c) < 0. Maximizing over c G C 
yields the bound Co<t>(x) + U*(<j>'(x)) < 0 and, therefore, W is a supersolution of 
the HJB equation. □ 

Lemma 10.3 If I25\) holds then W > 0 is a viscosity subsolution of S 1 3\) . 

Proof. Let x € int K and let tj> € C 1 (K ) DC 2 (x) be a function such that <j>{x) = W(x) 
and W < cf> on K. Suppose that the subsolution inequality for cf fails at x. Thus, 
there exists e > 0 such that C<f> < — e on some ball O r {x) C int K. By virtue of 
Lemma [8.31 (applied to the function <f>) there are to > 0 and rj > 0 such that on the 
interval ]0,to] for any strategy n e 

E (jf AT ~ + e~P T *cf>(y?£«)I {t Ar-<e } ) < 4>{x) - Tit, 

where r’ 1 ’ is the exit time of the process V x,7r from the ball O r {x). Fix arbitrary 
t G]0,to]- By the second claim of Lemma [9.11 there exists n G Aa such that 

W(x) < E (j? AT + e-P T W(V?£)I {tAT<9} ) + \r)t, 

for every stopping time r, in particular for tT 

Using the inequality W < <j> and applying Lemma f8.3l we obtain from the above 
relations that W(x ) < tf>(x) — (1/2 )r/t. This is a contradiction because at the point 
x the values of W and <j> are the same. □ 


11 Uniqueness Theorem 

Before formulating the uniqueness theorem we recall the Ishii lemma. 

Lemma 11.1 Let v and v be two continuous functions on an open subset O C R fi . 
Consider the function A(x, y) := v(x ) — v(y ) — \x — y \ 2 with n > 0. Suppose that A 
attains a local maximum at (x, y). Then there are symmetric matrices X and Y such 
that 

( n(x-y),X) G J + v(x), (n(x-y),Y ) G J~v{y), 


X 

0 



< 3n 



and 


(27) 
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In this statement I is the identity matrix and J + v and J~v are values of the 
set-valued mappings whose graphs are closures of graphs of the set-value mappings 
J + v and J— v, respectively. 

The inequality (1271) means that for any vectors x and y from R <; 

(Xx,y) - (Yy,y) < 3n\x - y\ 2 . (28) 

Of course, if v is smooth, the claim follows directly from the necessary condi¬ 
tions of a local maximum (with X = v"(x), Y = v"(y) and the constant 1 instead 
of 3 in inequality 071) ). 

The inequality (1271) implies the bound 

tr (A(x)X - A(y)Y) < 3n|,4 1/2 | 2 |a: - y\ 2 (29) 

which will be used in the sequel (for the proof see, e.g., Section 4.2 in {19]'). 

The following concept plays a crucial role in the proof of the purely analytic 
result on the uniqueness of the viscosity solution which we establish by a classical 
method of doubling variables using the Ishii lemma. 

Definition. We say that a positive function £ £ C\(K) nC 2 (intJf) is the Lyapunov 
function if the following properties are satisfied: 

1) £'{x) G int K* and Cq£{x) < 0 for all x G int K, 

2) £{x) —> oo as |a;| —> oo. 

In other words, £ is a classical strict supersolution of the truncated equation 
(without the term U*), continuous up to the boundary, and increasing to infinity 
at infinity. 

Theorem 11.2 Suppose that there exists a Lyapunov function £ and the Levy measure 
II is such that 

II (z : x + D £ z G dK}) =0 Vi G int K. 

Then the Dirichlet problem w, 0 has at most one viscosity solution in the class 
of continuous functions satisfying the growth condition 

W (*)/£(*)—>■ 0, |*| —> oo. (30) 

Proof. Let W and W be two viscosity solutions of m coinciding on the boundary 
dK. Suppose that W(z) > W{z) for some z £ K. Take e > 0 such that 

W(z) - W{z) - 2 e£(z) > 0. 

We introduce a family of continuous functions A n : K x K —> R by putting 

A n (x,y) := W{x) -W(y)~ ^n\x - y\ 2 - e[£(x) + £{y)\, n > 0. 

Note that A n {x,x) = Aq(x,x) for all * G K and zlo(*,*) < 0 when x G dK. From 
the assumption that the function £ has a higher growth rate than W we deduce 
that A n (x,y) —> —oo as |*| + |j/| —> oo. It follows that the level sets { A n > a} are 
compacts and the function A n attains its maximum on a compact subset of K x K 
which does not depend on n. That is, there exists ( x n ,Vn) G K x K such that 

An(x n ,yn) = An \= sup A„(x,y) > A := sup A 0 {x,x) > 0. 

( x,y)£KxK x£:K 
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All ( x n ,yn ) belong to the compact set {{x,y) : Ao{x,y) > 0}. It follows that the 
sequence n|a;n — y n | 2 is bounded. We continue to argue (without introducing new 
notations) with a subsequence along which ( x n ,yn ) converge to some limit (x,x). 
Necessarily, n\x n — yn\ 2 —> 0 (otherwise we would have Aq(x, x) > A). It is easily 
seen that A n —> Ao{x,x) = A. Thus, x is an interior point of K and so are x n 
and y n for sufficiently large n. Let j > 1 be the number of nonzero components 
x. Without loss of generality we assume that x J+1 ,... ,xd = 0 and for sufficiently 
large n the first j components of x n and y n are strictly positive. 

By virtue ofthe Ishii lemma applied to the functions v := IT —el and v : = W+e£ 
at the point ( x n ,yn ) there exist matrices X n and Y n satisfying ( 1271 ) such that 

(n(x n — y n ),X n ) G J + v(x n ), {n(x n - y n ),Y n ) G J~v(y n ). (31) 
Suppose for a moment that 

{n(x n — y n ),X n ) G J + v(x n ), {n(x n - y n ), Y n ) G J~v{y n )- (32) 

Using the notations p n ■= n(x n - yn) + e£'(x n ), q n ■= n(x n — y n ) — e£'{yn) and 
putting X n := X n +e£"(x n ), Y n := Y n —el''{y n ), we may rewrite the last relations 
in the following equivalent form: 

{p n ,X n )eJ + W(xn), (q n ,Y n ) € J~ W (j/„). (33) 

Since W is a viscosity subsolution, by virtue of Lemma[T2]there exists a function 
fn € Ci (A') n C 2 (x n ) such that f' n {x n ) = Pn, fn(xn) = X n , fn(x n ) = W(x n ), and 
W < f n < W + 1/n on K. Since IT is a viscosity supersolution we conclude in the 
same way that there exists a function f n € Ci(A') n C 2 (y n ) such that fh{yn) = qn, 
fn(yn) = Yn, fn(yn) = W (y n ) , and W - 1/n < fn < W on K. To deal with the 
nonlocal integral operator we take f n and f n having the structure given in ( 1181 ) 
with an appropriate choice of the cylindrical cutoff functions. We discuss details 
of this choice later. 

By definitions of sub- and supersolutions we have that 

F(X n , p n , T( f n , X n ), W(x n ), Xn) > 0 > F(\ n , qn,3-(fn, yn), W{yn), yn)- 

The second inequality implies that mq n < 0 for each m £ G = (—A) n9Ci(0). But 
for the Lyapunov function £'{x) G int K* when x G int K and, therefore, 

mpn = mq n + em(f'(x n ) + i'{y n )) < 0. 

Since G is a compact, Eq{p h ) < 0. It follows that 

F 0 {Xn,Pn,T{fn,Xn),W{x n ),Xn) + U* (p n ) > 0, 
Fo{Yn,q n ,F{fn,yn),W(yn),yn) + u*{q n ) < 0 . 

Recall that U* is decreasing with respect to the partial ordering generated by C* 
hence also by A'*. Thus, U*{p n ) < U*{q n ) and we obtain the inequality 


bn - — Fq {Xn, Pn, '-£{fn, Xn), IT {xn), Xn) Fq {Yn, qn, F{ fn, Pn) , IT {yn), Pn) ^ 0. 
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Clearly, 


<2 d 

bn = 2 ( al3x n x3 n X i 3 - yWnYij) + y (x l n - y n ) 2 

i,J = l 2=1 

-I/3n|a;n - 2 m| 2 - /3A n (x n , y n ) + X(f n - e£,x n ) - X(f n + e£,y n ) 

+e(C a£(x n ) + Co £{yn))- 

By virtue of (1291) the first term in the right-hand side is dominated by a constant 
multiplied by n\x n — 2/n| 2 ; a similar bound for the second sum is obvious; the last 
term is negative according to the definition of the Lyapunov function. To complete 
the proof, it is sufficient to show that 

lim sup (l(fn - e£, x n ) - X(f n + e£, %in)) < 0. (34) 

n 

Indeed, with this we have that lim sup b n < —/3A < 0, i.e. a contradiction arising 
from the assumption W(z) > W(z). 

In general, we cannot guarantee that Ol holds for arbitrary test functions f n 
and f n - That is why we choose them in accordance with the expressions given by 
Lemma 1151 with a = 1/n, i.e. 

fn(x n +h) := /S(*n + h)(a„,< (h) + W(x n + h) (l - (h )), Xn+hdK, 

fn(yn + h) := f°(y n + h)g nfi , a (h) + W{y n + h)( 1 - y n + h e K, 

where the linear space L := Im D% = {x € R d : x^ +1 = 0,..., x d = 0}, 

f°(x n + h) := [(W(x n ) + Q Pn ,x n (h) + r n {h )) V W n (x n + h)] A [W(x n + h) + 1/n], 
f°{y n + h) := \{W(yn) + Qq n ,Y n (h)-r n (h))AWn(yn + h)] V [W{y n + h) - 1/n]. 

Let 8 := (l/2)_R/(|a;| + R ) where R = d(x,dK) is the distance of the point x 
from the boundary dK. Then I(z,x n ) = 1 and I(z,y n ) = 1 for z € 0^(0) when n 
is sufficiently large. Indeed, when \x n — x\ < R/2 and \y n — 3i| < R/2 we have that 

|Xn + D Xn z - x\ < R/2 + \x n \\z\ < R/2 + (|J| + R/2)\z\ < R 

when |s| < <5, and the similar estimate holds for y n - 
We have: 


X(f„-e£,x n )-T(fn + e£,yn)= [ H n (z)n(dz)+ I H n {z)IT{dz) 

J{\z\<6} H\z\>S} 

where H n (z) := F n (z) — F n (z) with 

F n (z) := (fn - s£)(x n + D Xn z)I(z,x n ) - (w - e£)(x n ) - (fn - e£')(x n )D Xn Z 
= (f n - s£)(x n + D Xn z)I(z,Xn ) - (W - s£)( Xn) ~ n(x n ~ Vn)D Xn Z, 
F n (z) := (f n + e£)(y n + D Vn z)I(z,y n ) - (W + e£)(y n ) - (fn + e£')(y n )Dy n z 
= (fn + e£)(y n + D Vn z)I(z, y n ) - (W + e£)(y n ) - n(x n - Vn)D yn 


z. 
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Define also the functions 


G n (z) := (W - e£)(x n + D Xn z)I(z,x n ) - (W - e£){x n ) - n{x n - y n )D Xri z , 
G n {z) := (W + s£)(y n + D Vn z)I(z,y n ) — (W + e£)(y n ) - n(x n - y n )D Vn z. 


Put x' n := PlXu, y' n '■= PlVu- Considering only sufficiently large n, we may 
assume without loss of generality that (1/2)2® < x l n , yh < x 1 + 1 for non-zero coor¬ 
dinates of x and, therefore, the norms of the restrictions of the diagonal operators 
D x i , D y i on L and their inverses are bounded (even uniformly in n). Therefore, 
we may apply Lemma 17.21 if x has no zero components or its extension given by 
the accompanying remark and argue further supposing that 


Za n ,a' n (D Xn z) = Za n ,a'J D X ' n z) = tL,<( D v' n z ) = Vz G ^ 


According to our choice of 5 for z £ £> 5 ( 0 ) we have I(z,x n ) = I(z,y n ) = 1 for 
sufficiently large n and, as a consequence, the following easily verified identity: 


Gn(z)-G n (z) 


A n (x n + D Xn z, y n + D Vn z) - A n (x n , yn) + ^n\(D Xn 


Dy n )z\ 


2 


Recalling that the function A n (x,y ) attains its maximum at ( x n ,yn ) we get from 
here the bound 

Gn(z) - Gn(z) < i n\x n - y n \ 2 \z\~. 

Thus, on the set {z : \D x ^z\ > a' n } n 0^(0) we have that 


H„{z) = G n {z) - G n (z ) < in|x n - y n \ 2 \z\ 2 . 
On the set {z : \D x ^z\ < a n } D 0 , 5 ( 0 ) 


(35) 


Fn{z) = i(A - n D Xn z,D Xn z)+r n {D Xn z) 

+£^(ln) +£'(x n )D Xn Z+ -(£"(x n )D Xn z,D Xn z) — £{x n + -D x „z)J, 
Fn(z) = i( Y n Dy n z,Dy n z ) + r n (D Vn z) 

+£^£{y n ) + f!(y n )D yn z + i(/' (y n )D yn z, D Vn z) - £(y n + D Vn z) j. 
By the Ishii lemma 

(. X n D Xn z , D Xn z) - ( Y n D yn z , D Vn z) < 3n\D Xn z - D yn z\~ < 3n|x„ - y n \ 2 \z\ 2 . 


We take a n small enough to ensure that 

/ (r(D Xn z) + r(D Xn z)n(dz) < 1/n. 

■’ {\ D x' n z \< a ™} 

On the set {z : \D x i z\ > a n } n 0<s(0) the sequence of functions \H n \ is bounded 
by a constant. We choose a' n sufficiently close to a n to guarantee that 
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The expressions in parentheses in the above formulae for F n (z) and F n (z) (residual 
terms in the Taylor formula for the smooth function £) are bounded by a constant 
times \z\ 2 and converge to the same limit as n —> oo. 

Summarizing the above facts we conclude that 

limsup / H n (z)II(dz) = 0. 
n -'{|z|<<5} 

Since the continuous function W and l are of sublinear growth and the se¬ 
quences x n and n(x n — yn) are converging (hence bounded), the absolute value of 
F n is dominated by a function c(l + \z\). The arguments for — F n (z ) are similar. 
So, the function H n is dominated by a function of sublinear growth. 

Put Z := {z : x + D.^z G dK}. By assumption, Z is 17-null. li z ^ Z and n is 
large enough, we have I(x n ,z) = I{y n ,z) = I(x,z). It follows that 

H n {z) = ( Gn(z ) - Gn(z))I(x, z) - n(x n - yn){D Xn - Dy n )z( 1 - I(x,z)) 

~{W - fn) {Xn + D Xn z)I(x, z) + (W — fn){yn + Dy n z)l(x,z) 

+ (W + e£)(y n ) - (W -e£)(x n )) (1 -I(x,z)) 

< ( A n (x n + D Xn z, y n + D Vn z) - A n (x n ,y n )) I(x, z) 

+ n\x n -y n + (D Xn - D Vn )z\ 2 - ^n\x n - yn\ 2 ^j I(x,z) 

- n(Xn - yn){D Xn - D y „)z -I- 2/n 

F(W + e£){y n ) - (W - e£)(x n ))(l ~ I(x, z)). 

The first term in the right-hand side is negative as A n (x,y) attains its maxi¬ 
mum at (x n ,yn)- We only need to consider the second term when I(x,z) 0. In 
this case, we may combine it with the third term and conclude that the sum is less 
than (l/2)n|a;n — yn | 2 —► 0. The last term converges to zero because of continuity. 
Thus, 

limsup / H n (z)II(dz) = 0. 

n -l{\z\>5} 

However, our reasoning is based on the assumption (El while we know only 
El- Fortunately, using the definitions of J + v{x n ) and J v(y n ) we can replace the 
objects x n , yn, X n , Y n by their approximations x n , yn, X n , Y n approaching rapidly 
the initial ones and for those El hold. Repeating the arguments and controlling 
the approximation errors, we get the same contradiction. □ 

Remark 1. Note that the definition of the Lyapunov function does not depend 
on U and hence the uniqueness holds for any utility function U for which U* is 
decreasing with respect to the partial ordering induced by K*. However, to apply 
the uniqueness theorem one needs to determine the growth rate of W and provide 
a Lyapunov function with a faster growth. 


12 Existence of Lyapunov Functions and Classical Supersolutions 

In this section we extend results of [18] on the existence of Lyapunov functions 
and classical supersolutions to the considered case of nonlocal operators. 

Construction of Lyapunov functions. 
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Let u G C(R+) n C 2 (R+ \ {0}) be an increasing strictly concave function with 
w(0) = 0 and u{ oo) = oo. Introduce the function R := —u' 2 /{u"u). Assume that 
R := sup z>0 R{z) < oo. 

For p G K* we define on K the positive function / G Ci(K) n C 2 (int K by 
putting f(x) = f p (x ) := u(px). If y G AT, then yf(x) = ( py)u'(px ) > 0. 

If p G int A'*, then for any x,y £ A'\{0} we have the strict inequality yf'(x) > 0 
implying that f'(x) G int A"*. Thus, for p G int A'* the function / is a Lyapunov 
function provided that the inequality Cof(x) < 0 is satisfied. We show that under 
some mild conditions this inequality holds for sufficiently large ft. 

Put k p := 0, r/p := 0, if II = 0 and 


u (px) 

k p := sup ——r- 

aiGint K UypX) 


\p\\x\, r/p := Kp / \z\n(dz) 

J{\z\> K ^} 


otherwise. Define also 


r/p := - sup 


(/i(z),p ) 2 


2 xGi ItK (A(x)p,p) 1 ^ A ^^- 


Recall that A{x) is the matrix with A lJ (x) = a lJ x l x J and the vector y[x) has 
the components p l x l . 

Note that if k p < oo, then r/ p < oo (as we assume that f \z\ 2 A \z\II(dz) < oo). 

Example. Let u(z) := z p /p where p g] 0, 1[. Then R = R(z) = p/( 1 — p) and, for 
p G int A'* 

. |p||x| 

Kp < p SUp - < OO 

xGif\{0} P X 

(strictly positive function y i-a py on the compact K C\ {y: \y\ = 1} attains its 
minimum). 

Proposition 12.1 Let p G int A'. If k p < oo and /3 > fj p R + r/ p T max,; \pi\Kp, then 
f p is a Lyapunov function. 

Proof. Let x G int K. Recall that 

£(/,*) := j [(/(* + D x z)I intK {x + D x z) - f(x)) - D x zf'(x)]n(dz). 

If x + D x z G int A', then the integrand defining T(f, x) has three nontrivial terms 
and we have by the Taylor formula (in which i9 G [0,1]) that 

(f(x + D x z) - f(x) - D x zf(x )) = i u”(px + DpD x z)(pD x z) 2 < 0. 

If x + D x z (f int K. then the integrand is reduced to two terms. Moreover, for 
\z\ < 1/k p we have the bound 

\D x zpu (px))\ < |«||p||a;|M (pa;) < u{px) 


implying that 


/(*) - D x zf'{x)) = - u{px ) - D x zpu (px)) < 0. 
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We obtain from here, taking into account that u(px) > 0, the bound 

Af,x) < u{px)\p\\x\ I hnt k(x + Dxz)n(dz) < Tjpu(px). 

J{\z\>1/k p } 

Suppose that (A(x)p,p) ^ 0. Isolating the full square we obtain that 
£ o/(*) = \ (A(x)p,p)u"(px) + 2{p(x),p)u\px) + ^fypl) l"(px) 
1 (^(x),P) 2 


+ 2 (A(x)p,p) 
1 (p(x),p) 2 


R(px)u(px) + T(f, x) — /3u(px) 


- 0 , ., s . R(px)u(px) + v P u(px) - /3u(px ). 

2 (A(x)p,p) 

It follows that Cof{x) < 0 if /3 > fj p R + p p . 

Of course, if (A(x)p,p) = 0 we cannot argue as above. In this case 

Cof(x) < (p(x),p)u (px) +p p u(px) - /3u(px). 


Taking into account that 


(p(x),p)u'{px) 

blip . \ 

xGint K U[px) 


^ I I u'ipx) | 11 . 

< max \pi\ sup ——dp i 

-- n, ,iVX' 


r-Ci nt IT U 


max |pi | ftp, 
i 


we get that £o f(x) < 0 , if 

ft > Vp + max \pi\Kp. 
i 

Combining these two cases, we get the result. □ 

Remark. An inspection of the above arguments shows that one can get that the 
f p is a Lyapunov function for 


> sup 
fcGint K 


f 1 (u(x),p) 2 , . 

\ 2 (A(x)p,p) R ^ I {{A{x)p,p) 1 lo} + 


(p(x),p)u'(px) 
u(px ) 


\{A(x)p,p)A0} \+Vp- 


Of course, such a bound is less tractable than that given above. 

Construction of classical supersolutions. 

Similar argments are useful in the search of classical supersolutions for the 
operator £. Since £/ = £o/ + U*(f'), it is natural to choose u related to U. 
For a particular case, where C = R+ and U(c) = u(eic), with u satisfying the 
postulated properties (except, maybe, unboundedness) and assuming, moreover, 
that the inequality 

u* (au (z)) < g(a)u(z) (36) 

holds, we get, using the homogeneity of Co, the following result. 

Proposition 12.2 Let p € intff. Suppose that IS6\) holds for every a, z > 0 with 
g(a) = o(a ) as a —t oo. If k p < oo and ft < rj p R + g p + max^ \[m\k p , then there exists 
a o such that for every a > a o the function af p is a classical strict supersolution of 

m 

For the power utility function u(z ) = z 1 / 7 , 7 €]0,1[, we have: 

u*(au(z)) = (1 — 7 )a y ^ 7 ~ 1 ' > u(z). 


Therefore, the inequality (1361) holds with g(a ) = o(a), a —> 0. 
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